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FOREWORD 


This is the final report submitted on Contract NAS5-2108 
entitled, "Hybrid (Hydrazine) Resistojet Development " The program 
originated in the Auxiliary Propulsion Branch of the NASA/Goddard 
Space Flight Center. Mr. Robert A. Callens was the Technical 
Officer for the NASA/ Goddard Space Flight Center. The Project 
Director at Avco Corporation was Mr. T. K. Pugmire. Other 
participants in the program were W. Davis, T. O'Connor, R. 
and R Shaw 


Olson, 



I Summary 


This report discusses work performed under NASA Contract 
NAS5 -21080, "Hybrid (Hydrazine) Resistojet Development " The ex- 
perimental evaluation of several hydrazine thrusters which were ignited 
by electrical heaters within the thrust chamber indicated that high specific 
impulse (' v 240 sec.) could be achieved during steady state operation with- 
out power addition, except for that required for ignition purposes Restart 
capability was demonstrated with the various thrusters investigated 

Pulsed operation of preprototype thrusters indicated that the 
specific impulse obtained for pulses on the order of 1 sec duration de- 
creased much more rapidly than that observed with nominal 5 lbf catalytic 
hydrazine thrusters. A limited amount of experimental data for pulse 
widths of less than 0.1 sec , however, demonstrated that performance 
comparable to that of catalytic engines can be achieved at duty cycles 
greater than one percent. 
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II 


Conclusions 


The experimental data gathered during this study indicate the 
feasibility of developing a hydrazine thruster in which propellant 
decomposition is initiated by thermal energy stored within the thruster 
components or by an electrical heater which would be employed in the 
initial start-up of the unit Specific impulses well m excess of 200 sec 
can easily be obtained over a wide range of thrust levels Experimental 
evidence demonstrates that this type of performance can be achieved 
even at a thrust level of 5-10 millipounds simply by reduction of heat 
losses from the thruster 

Operation of a hydrazine resistojet thruster at specific impulses 
on the order of 75 percent of steady state values can be achieved for 
short pulses The length of the pulse for which this type of performance 
can be achieved is governed by the ease with which energy can be trans- 
ferred from thruster components to the propellant 
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Ill 


Analysis 


Hydrazine is an unstable chemical compound as is demonstrated 
by its standard free energy (AF°£2gg=38. 02 kcal/mole){l) which makes 
it suitable as a monopropellant The decomposition of hydrazine may 
be written as 

M a H 4 - ~ NH -3 + ± -AH=33 8 k cal/ mole (1) 

or 

N 2 w 4 2 ■** N j> -Ah= 22 75 k cal/mole (2) 

In actual practice, however, the decomposition occurs m accordance with 
stoichiometric equation 

n & h 4 * nh 3 + axu^ > (i+2X)n 2 (3 ) 

3 3 

where the quantity, x, reflects the amount of deviation from the most 
energetic chemical reaction (eq 1) Experimentally, it has been found 
that this fraction, x, is a function of the temperature at which the re- 
action occurs (2). 


The actual decomposition of hydrazine occurs by a series of steps 
which are considerably more complex than the overall stoichiometric 
relationships (eqs 1-3). For example, "Eberstem and Glassman^), m 
evaluating their experimental homogeneous reaction rate data gathered 
over the temperature range of 750 to 1000°K, proposed the following de- 
composition . 


Initiation Reaction. 

N 2 H 4 +M Z M 

Propagation Reactions 

N*W 4 + NH, NH 3 + M 2 VA s 

+ ^ N a + H z 4- H + M 

NHj+NHj 

Branching Reactions 

+ M ► NH T N-H a 4- M 

+ NH N 4- 


K= > Q 1 ^ exp feO,OQO 

' RT 


K = 

io® 

exp(- 

TOOO\ 

(5) 



RT ) 

K = 

VO* 

exp{- 

20,000\ 
RT J 

(6) 

K ~ 

Kf 

exp(- 

■7O0O\ 
RT J 

(7) 

K = 

\o a 

9 exp( 

- t8,OOQ''y 
RT / 

(8) 

K = 

io' 4 exp ^ 

-io,ooo\ 
RT ) 

(9) 
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Termination Reactions 


NiH 5 NiH 3 +N z + H t 

N z H 3 + H 3 > 2NH 3 +Na. 

N x H 3 *** ^ + 2. H 

NHa-v-KHa ^ N 7 .H 4 


»z s 

K- \0 

(10) 

X 

n 

0 

N 

(11) 

IS 

K- \0 

(12) 

K- 'O 

(13) 


At low temperatures, the predominant propagation reactions are 
those described by equations 5 and 7 making the final products formed 
in the termination reactions mostly ammonia and nitrogen However, 
at higher temperatures, the propagation reaction which produces nitrogen 
and hydrogen {eq 6) and the associated termination step (eq 12) become 
more important leading to larger and larger amounts of hydrogen m 
the final products. Notice that the reaction scheme above does not in- 
volve the homogeneous decomposition of ammonia 


nh 5 


J_ Si, •+ 3 H, 
2 z a • 


{14) 


which has been found to be slow(^) 

= - 10 s 64 exp [ m] 

which implies that in a volume typical of the thrusters considered m 
this program that a residence time of 3 x 10 "^ sec at a temperature of 
3000°R would result in only a 5 percent ammonia decomposition. 


man^) 


The overall kinetics which were obtained by Eberstem and Glass- 

± [N, n»l = - K [N 4 H 4 ] [m] 

dt < 16 ) 


is second order due to the low hydrazine pressures at which these ex- 
periments were performed. However, at high pressures the decomposition 
of hydrazine is found to be a ummolecular reaction^ ' expressed as 

L^U*] = -K H 4 ^] (17) 

dt 

with the transition occurring at pressures m the vicinity of one atmosphere. 
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Since a hydrazine thruster would operate with chamber pressures 
m excess of atmospheric, the kinetics of hydrazine decomposition are 
taken to be first order. The amount of ammonia present in the products of 
the homogeneous decomposition must, however, be taken to be a function 
of temperature as shown in Fig. 1 which includes the experimental data of 
Gray et al, McHale et al^) and Michel and Wagner^) It can easily 
be seen that the amount of ammonia present m hydrazine decomposition 
products corresponds to x< 0 5 m the stoichiometric equation 

n*h 4 -|(i-x)nh 3 + 2xh 1 + ± (i + ax^Nj. (18) 

for all temperatures less than 2500°R which correspond to the adiabatic 
flame temperature at x = 0 5 


As has been cited previously, the further decomposition of ammonia 
by dissociation is a slow homogeneous reaction and would not be expected 
to alter the composition of hydrazine decomposition products. It is well 
known, however, that many metals act as catalysts for the heterogeneous 
decomposition of ammonia (^) The list of catalytic metallic surfaces 
includes iron, molybdenum, tungsten, rhenium, ruthenium, osmium, 
copper, cobalt, mckel, rhodium and platinum In general, the decomposition 
rate per unit area of exposed surface is given by 

T, = ke *p[-| F )[NH,] X [H z T' (WJ 


with the exponents x and y varying with material The values of the 
exponents typically vary from both x and y being zero for tungsten and 
molybdenum to x = 0. 7 and y = 1 4 for rhenium. Typical activation energies 
(E) required for the process are generally m the range of 40+5 k cal/mole 


and the frequency factor, A, is on the order of 10^ 


Comparison of the 


catalytic decomposition rate (equation 10) with the homogeneous reaction 
rate (equation 15) yields 


r c _ 


exp (-4 Z,soq/rt ) 

exp 7^500/RT ) 


= \o 


exp ( 37,000/RT) 


which clearly demonstrates that, at the temperatures expected within the 
hydrazine thruster, the heterogeneous decomposition may be considerably 
more important than homogeneous gas phase reactions. 


The amount of hydrazine which can be reacted within a given volume 
may be easily obtained from data reported m the literature and analysis 
Kinetic data at pressures in excess of atmospheric pressure and at 
temperatures greater than 1000°K are illustrated m Fig 2 The data ob- 
tained by McHale et al^ are applicable over the temperature range of 975 
to 1540°K while the curve obtained from Eberstem and Glassman^) is a 
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fig 2 WYDRAIINE DECOM POSITION KINETICS 
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correlation of original data and results published m the literature. 
Since the latter curve can be applied over a wider temperature range, 
it will be employed m the subsequent discussions 


A material balance for hydrazine contained in any volume may be 
written as 


■vin ( W L - w 0 ) = - rV 


where it is assumed that the fluid within the volume is completely mixed 
witn a mass fraction equal to that of the fluid leaving the volume through 
the exit nozzle The reaction rate of the propellant is given by first 
order kinetics as 


Y* = - kpw 0 e/rT^ 


with being the average concentration within the system. Combining 
these two relationships, we get 

«*p(-e/rt) 

with the average density obtained from the equation of state 

p •= PM 

■RT 

The mass flow leaving the system is simply 


PA., f 2 \*< 

J RT ' H+l / 


Hence, 


Wi i — __ 

A, 


y K exp(-g/RT) 

- w W)W 


Since by definition 

L*= V. 
A» 


* • 




L* K exp(- E/RT) 
■ " ' " 5 SL 

1 
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Notice that the quantity lT / a is a characteristic "residence time" m 
the system Equation 30 may be rearranged to give the following 
relationship for the ratio of outlet ^.nd inlet propellant concentrations 










(31) 


From the first order rate expression, we may define a half life for the 
leactant at any average temperature since 


^ w - W exp f T.^L-N 

at ' 

or upon integrating m the case of a constant temperature system 

ln ( >0 = - K [*^(-^r)] t 

The time at which the concentration decreased to one -half of its initial 
value is then 

X _ 0.693 

vw ~ - ■ ■ . ■ ■-> 


(32) 


(33) 






(34) 


which may be substituted into equation 31 to give 


J±L 


Wo _ 






(SJ\ aot-o Jr O 693 ( -^=- ) 
‘ i / ' CXty 


V'frH ! v atyj) 

The quantity (IT / a t 1/2) represents a dimensionless time for the system 
Large dimensionless times imply a large fraction of propellant is de- 
composed. This can best be accomplished by increasing the temperature 
of the system since 


a ■fc'A 


oC \L ex P (-JL\ 


(35) 


with the exponential increasing much more rapidly than T 


V2 


The ratio of outlet to inlet mass fractions predicted by equation 31 
are presented in Fig 3 Over the range of mass fraction ratios less than 
0.25 the conversion is well represented by 

o^5 


vy. f kiJ* ^cp /- e \\ - 

Kia 1 


0 493 


(36) 


Remembering that 




'I 


M 


’/Z. 


(37) 
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and if it is assumed that the average molecular weight of the combustion 
products is 15.2 (corresponding to 95 percent conversion), the following 
relationship is obtained 


We 

W v 



(38) 


At a characteristic length of 1 ft. and a temperature of 2000°R, it is 
found from Fig . 4 . that 

K <5xp ( - H \ = 158X10 3 

V KT ) 



(39) 


(40) 


Thereby indicating that the ratio of outlet to inlet mass fractions for this 
particular case is 0 242. Since the reaction rate increases exponentially 
with temperature, increasing the temperature within the volume to 2500°R 
decreases the mass fraction ratio to a value of 4 2 x 10“ ^ for L* equal to 
1 ft. Increasing the characteristic length to 10 would further reduce the 
mass fraction ratio, i e , increase the conversion of propellant to reaction 
products as shown in Fig 4. 


Since the concept of a thermally ignited hydrazine resistojet 
also has the capability of using power addition to enhance the performance 
of the system, analyses were conducted to evaluate the magnitude of the 
increase which could be obtained by means of simple energy balance 
relationships. 


The chemical power input to such a thruster is given by 

£ 


^ + jCp dT j 


(41) 


For adiabatic operation of the thruster, the power in the exhaust 
gases is equal to the chemical power input 

^chem ® W J Cp t c ^ 

l ° 

which allows the adiabatic flame temperature to be established for any 
inlet temperature T^ since 


(42) 


I W 'I ^ 


dT = aH. 




I 


(43) 


for any set of reaction products 
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W c /W t , UN REACTED FRACTION 


PHG 4- HYDRAZINE DECOMPO£\T\OH DEPENDENCE 
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The decomposition reaction of hydrazine is given by 

*. 4 (\-x) N W* + 2X H z + ( i +-ZX) N, (44) 

i 3 3 3 3 3 

where X is the fraction of ammonia dissociated. The energy released m 
this exothermic reaction of liquid hydrazine to gaseous products^) is 

AH- = - 1505 +• 82.S X BTU / LB (45) 

This is equivalent to the energy release of 33 k cal/g-mole reported m 
the literature(9) for the gas phase reaction with 25 percent of the ammonia 
dissociated. 

* 

Also, the adiabatic flame temperature at various ammonia dissociation 
fractions was determined by using thermodynamic data for the various chemical 
species as listed in the JANAF tables ^ and are illustrated m Fig 5 m the 
form of power addition per unit mass flow (watts-sec /,/flb m ) as a function 
of temperature with dissociation fraction as a parameter The adiabatic flame 
temperature can be determined directly from Fig. 5 since this case corres- 
ponds zero power addition The amount of ammonia dissociation as re- 
ported m the literature^’ ^ ’ 6>is also included m Fig. 5 In the case where 
the reaction zone is adiabatic and the decomposition occurs homogeneously, 
these data indicate that the stoichiometry of the reaction is 

OT7NH, +0 84 +061 (46) 

* * *3 3 3 

and the adiabatic flame temperature is 2400°R In practice, however, 
the amount of dissociation will be less than that indicated by equation 46 
(X=0 42) since there will be some energy losses from any actual system, 
and heterogeneous reactions are expected to play some part in the de- 
composition 

The performance enhancement which can be obtained can be 
determined from the variation m specific impulse with temperature The 
specific impulse of a hydrazine thruster was determined analytically from 
the relationship 
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POWER ADDITION CLOGS), WATT/C IB m /SEC> 


PIG 5 EFFECTS OFPONWER At)WT\0^ ON TEMPERATURE AHt> 





with the molecular weight (M) being a function of the ammonia dissociation 
fraction and the specific heat ratio as a function of the amount of ammonia 
present m the decomposition products as presented in Ref 10 The 
specific impulse for a thruster having an area ratio of 80 1 is shown in 
Fig 6 for various gas temperatures. The theoretical performance 
curves demonstrate that for a completely adiabatic system the specific 
impulse would be m the rate of 250 sec. for x = 0 42 {homogeneous 
hydrazine decomposition) to 262 sec for x = 0. 0 (heterogeneous hydrazine 
decomposition) . 

The performance curves can be used to demonstrate the importance 
of further ammonia decomposition on the efficiency of power addition in 
increasing specific impulse of a thruster. For example, consider a unit 
where the decomposition of hydrazine occurred m an adiabatic homogeneous 
manner (X=0 42) Following the decomposition the thermal losses were 
such to produce a reduction of the gas temperature to 1800°F and, hence, 
yield a specific impulse of 220 seconds rather than the theoretical value of 
250 seconds If the impulse is to be increased by power addition to 260 
seconds and there is no further ammonia decomposition, the total power 
addition is 0 55 watts/ (>t<lb m /sec) However, if in the power addition pro- 
cess all the remaining ammonia decomposes the power requirement is 
increased to a theoretical value of 0-92 watts /(*lb m / sec). This simple 
exercise clearly demonstrates the advisability of limiting the amount of 
ammonia decomposition within the thruster if power addition is contemplat- 
ed. 


The increase in thrust which is achieved through power addition 
is shown m Fig. 7 The preparation of this illustration it was assumed 
that the thruster was capable of operation at a specific impulse of 180 sec 
and a thrust level of 10 millipounds with no power addition, and that the amount 
of ammonia in the exhaust was given by the ammonia "dissociation" fraction 
as a function of temperature A second power addition curve is also pre- 
sented which depicts the enhancement achieved with no change in composition. 
In the region where thrust increases linearly with power addition, a 1 milli- 
pound thrust enhancement requires 18 5 watts of power for variable com- 
position while the constant composition curve indicates that 14 5 watts of 
power are necessary to produce the same result 

These analyses clearly demonstrate the feasibility of the hybrid 
hydrazine resistojet concept Reasonable performance would be expected 
with a device which could operate at a specific impulse of 200 sec. at a thrust 
level of 20 millipounds This would be achieved with a thrust thermal 
efficiency on the order of 65 percent at a temperature m the vicinity of 
1700°R. Kinetics indicate that with a thruster m which propellant and re- 
action products were well mixed that 95 percent of the propellant would be 
reacted in a characteristic length of 10 ft. 
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IV 


Thruster Evaluation 


During the course of this study, several thruster configurations 
were employed to evaluate propellant injection schemes, the effects of 
power addition and to obtain thruster performance data. 

A) Laboratory Thrusters 

The devices employed in the earlier phases of the program were 
relatively simple units as shown schematically in Figs 8-10. Basically, 
these early thrusters consisted of a quartz tube mounted between a 
nozzle assembly and base plate assembly A heater coil of molybdenum 
or rhenium wire was employed as the means of raising the temperature 
of the thruster to levels sufficient to provide energy for propellant 
ignition. Several variations in the method of propellant injection were 
examined with this basic unit. These injection schemes are illustrated 
m Figs. 8-10. 

Injector configurations "A" consists of a simple porous plug 
mounted m the base plate with the thrust chamber being completely 
empty This configurations would be expected to have a minimum thermal 
inertia 


Configuration "B" thrusters were essentially identical to con- 
figuration "A" with the exception that porous ceramic rod (0 187 m dia 
x 1 0 m. long) was placed m the thrust chamber The heater coil m 
this configuration was in close contact with the rod of calcia stabilized 
zircoma 

In the configuration "C" thruster, the ceramic rod was mounted 
in the base plate and propellant flowed through the calcia stabilized 
zircoma by means of a 0. 030 in dia. passage located on the axis of the 
rod The ceramic rod had an external diameter of 0 187 in. and a length 
of 1. 00 in. The passage on the axis of the rod extended to approximately 
0 090 m. from the downstream end of the rod where it was sealed by 
ceramic cement. 

The last configuration (configuration "D") was nearly identical 
to "C" except that a 0.025 m o d. by 0 015 in. l d. tube was inserted 
into the ceramic rod for a distance of 0. 8 m from the upstream end ^here- 
by moving the point at which propellant enters the thrust chamber away from 
the thruster base plate. 

A limited amount of exploratory data was gathered with the configurations 
"A" and "B" while a larger effort was expended in evaluating configurations "C" 
and "D" Primary stability comparisons were made at zero power addition 
following the initiation of decomposition. The comparisons were made over 
a range of flow rates of 10“^ to 5xl0“^lb/sec (corresponding to a thrust level 
of 2 to 10 millipounds at a specific impulse of 200 sec ) In all cases, the 
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combination, of the orifice /porous injector (configuration "C") was 
superior to configurations "A" and "B" with respect to stable operation 
and decomposition efficiency as measured by characteristic velocity 
Samples of nozzle box and injection pressure responses for all con- 
figurations are shown in Figs. 11 , 12 and 13 at various mass flow rates 
at zero power input All configurations utilized the same exit nozzle 
and inlet fixture. The nozzle diameter was 0 017 in. which is a 
reasonable selection for a 10 millipound thruster In all cases, con- 
siderable conductive losses can be expected due to the intimate contact 
of the inlet and nozzle blocks with large heat sinks. These losses 
should be the same for all configurations so that good comparisons 
between competitive configurations may be made 

It is immediately obvious that the configuration "A" is 
extremely unstable at all conditions . The addition of the zirconia 
flame stabilizer greatly improved performance as can be seen by 
the pressure responses for configuration "B". When the performance 
of configuration "C" and "B" are compared, it is noted that the de- 
composition process for "C" is slightly more stable than "B" and at the 
highest flow rate, the characteristic velocity is 20% higher for the "C" 
configuration 

The large pressure fluctuations for the "A" configuration were 
primarily caused by flow variations through the porous injector. Once 
the thermal decomposition was initiated and the power turned off, the 
hydrazine was observed reacting on the exposed surface of the porous 
injector The reaction zone was a small area which was appreciably 
hotter than the mean temperature of the element. In addition to non- 
uniform surface temperature, the hot spot was observed to move around 
the injector, tracing a circle. This phenomenon was accompanied by 
rapid fluctuations m the temperature of the thruster heater wire pre- 
sumably caused by hydrazine passing through the porous element without 
reacting, and then decomposing on the wire. It can be assumed that some 
fuel was passed through the nozzle without decomposing at all as evidenced 
by the low mean pressures 

A reasonable explanation of the cause of these instabilities is as 
follows. When hydrazine vapor is entering the thruster at the burning 
velocity, we may write 

p = W (48) 

where the average density may be taken as that at the normal boiling 
point of the propellant (238°F). 

5 = PM (49) 

KT 
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For a burning velocity of 185 cm/secW, and a mass flow rate of 5x10-5 
lb/ sec, the minimum diameter of a cylindrical porous element is 


\f - 4wRT 

or 


(50) 


D=0 156 in. 


Notice that if the gas temperature is increased to 1250°F with the corres- 
ponding decrease m density, the minimum diameter is decreased to a 
value of 0. 024 m. 

Since a porous element similar to that used m configuration "A" 
and "B" thrusters cannot be expected to produce a uniform velocity 
stream, decomposition of the propellant would occur on those regions of 
the surface where the velocity is less than the burning velocity. These 
regions would tend to be localized hot spots due to heat transfer from the 
reaction zone to the porous element. However, as the temperature in- 
creases a critical value is reached when the local velocity is greater than 
the burning velocity and then the flame lifts off the surface. Furthermore, 
it is evident that an increase in the temperature of porous element would 
result in fluctuations m the mass flow through the system if the pressure 
drop across the element remains constant. This effect is caused strictly 
through the increase in viscosity with increasing temperature Hence, 
we may conclude that when the size of the porous element is close to 
that given by equation 50 than an unstable situation may result, e.g. the 
large scale pressure fluctuations observed with the configuration "A" 
prototype thruster. Smaller scale fluctuations are expected with con- 
figuration "B" since the ceramic stabilizer would tend to dampen the 
pressure fluctuations 

Stability and performance considerations indicated that configurations 
"C" and "D" (which was not evaluated m the experiments cited above) to 
be most satisfactory and further evaluation of these devices was merited 
Steady state operation of the configuration "C" thruster was examined over a 
flow rate range of 10“5 to 5 x 10 - ^ lb/sec; experiments were performed m 
both the powered and unpowered modes of operation These experiments 
were conducted with the thruster exhausting into vacuum with no provision 
for direct measurement of thrust, hence, the basic performance parameter 
employed m these evaluations was the characteristic velocity of the exhaust 
stream 


C* = 

rr> 


(51) 
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where the pressure, P]sj, is that measured in the thrust chamber, At 
is the nozzle throat area, and m is the mass flow through the system 

Propellant was supplied to the thruster by the hydrazine feed 
system shown schmetically m Fig. 14. The basic components of the 
feed system were a propellant tank having a capacity of 1000 ml which could be 
pressurized with inert gas to 180 psig, a flow metering system, a flow 
regulating valve, and a solenoid control valve. The flow metering device 
consisted of an orifice having a nominal diameter of 0 0043 m. which 
was instrumental with a differential pressure transducer (Dymsco Model 
DPT 69CB-10) which had a linear output over the range of 0-5 psid A 
schematic diagram of the propellant feed system is shownin Fig. 14 A 
series of calibration experiments were performed with the system by 
allowing hydrazine to flow through the orifice for periods ranging from 
180 to 240 sec. The total mass flow was collected in a beaker and weighed 
on an analytical balance having an accuracy of +0. 0005 grams. The re- 
sulting mass flow rate was correlated with the millivolt output of the 
pressure transducer as shown in Fig. 15 The data are well correlated 
by the equation 

>m = 3 75 x 10 5 ( E mv ) l© m /SBC * (52) 

with all data falling within 5% of this relationship 

Performance data were gathered with a configuration "C" thruster 
m both powered and unpowered operation over the flow rate range of 
1 5-5. 0xl0~5 lb m /sec. These data are summarized m Figs. 16-18 
in the form of characteristic velocity as a function of propellant flow rate 
with power applied to the heater coils as a parameter. Data obtained at 
a flow rate of 5 x 10~5lb m /sec are also presented m Fig 19 with characteristic 
velocity, specific impulse, and thrust as a function of power to the heater 
coils. The thrust was calculated from the experimental data using an 
assumed thrust coefficient of 1 65. It is seen from these results that a 50 
watt power input increases the thrust level from 7 2 x 10 - ^lb£ at zero 
power input to a level of 8 5 x 10“^lb£ This corresponds to a power re- 
quirement of 38.4 watts per millipound of thrust. When the .experimental 
value is compared with the theoretical thrust increase (Fig 7), it is found 
that the power addition process has an efficiency on the order of 50 percent 
when the variable composition thrust enhancement curve is used as a base 
A comparison of the calculated specific impulse with the theoretical pre- 
dictions at an assumed ammonia dissociation fraction corresponding to 
the exit gas temperature indicates a gas temperature at the nozzle of 
approximately 900°R whereas a completely adiabatic system would produce 
gas at 2400°R. The thermal efficiency of the system is also on the order 
of 50 percent. 
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A cursory examination of this system indicates that the thermal 
efficiency of the entire system would be expected to be quite low. For 
example, energy released by decomposition within the porous ceramic 
cylinder can be conducted into the base plate or may be lost by radiation 
from the external surface of the porous ceramic to the surroundings 
through the transparent quartz tube. Furthermore, considerable 
energy may be lost through heat transfer from the decomposition pro- 
ducts as they flow through the nozzle box assembly. It would be fortuitous 
if the system utilized in these preliminary experiments had a thermal 
efficiency much greater than 0 3-0.4. If a similar estimate is made of 
the efficiency of energy addition by means of the heater coils, it is 
immediately obvious that at least 50 percent of the power is lost by 
radiation to the surroundings. 

Visual observation of the thruster through the transparent 
quartz walls was conducted during all these experiments It was ob- 
served that at low propellant flow rates the decomposition of propellant 
occurred primarily m the region close to the thruster base plate. As 
the flow rate is increased the reaction zone became more extensive and 
extended further m the downstream direction. 

In spite of the deficiencies noted above, the configuration "C" 
thruster demonstrated restart capabilities without power application 
For example, with the thruster operating at steady state, the propellant 
valve was closed for periods ranging from 30 to 200 sec. Upon reopen- 
ing the valve, propellant ignition was found to occur The difficulty en- 
countered in the restart was dependent on the length of the off period and 
the magnitude of the propellant flow rate as is shown by the experimental 
records obtained m these investigations (Figs. 20-33) At low flow rates, 
rise times in the range of 5 sec. were achieved for off times on the order 
of 70 sec. 

In several instances, there was evidence of partial flooding of 
the thruster (Figs. 27, 28 ) Immediately upon opening the propellant 
valve, the nozzle pressure was found to rise sharply for a period less 
than 1 sec. For the next 10 to 15 sec. , the propellant flow rate was in 
excess of the final steady state flow rate, and there is an abrupt change m 
the slope of the pressure rise curve. After this period, steady state 
operation was achieved. 

Attempts were made to obtain similar data with this configuration 
thruster mounted entirely within the vacuum system as a means of in- 
creasing thruster efficiency by elimination of convective losses from the 
thrust chamber external surfaces The problems encountered were as 
follows 
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Fig. 21 Configuration "C" Laboratory Thruster Transients (20 watts) 
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Fig 22. Configuration "C" Laboratory Thruster Transients (No Power) 
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Fig. 23 Configuration "C” Laboratory Thruster Transients (No Power) 
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Fig. 24. Configuration "C" Laboratory Thruster Transients 
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Fig 31. Configuration. "C" Laboratory Thruster Transients 
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(i) Propellant leakage at the base of the porous zirconia 
injector, thereby allowing a porticn of the propellant to be 
introduced into the thruster on the external surfaces of the 
zirconia. This produced fluctuations in the pressure m the 
thruster nozzle. The area of this leakage is designated as 
area "A" m Fig. 34 of this report. 

(n) During operation temperatures, the thruster achieved levels 
sufficient to cause decomposition of the O-nng seals, particularly 
that seal at the nozzle end of the thrust chamber The area of 
this problem is designated as area "B" in Fig. 34. 

( 111 ) Increased temperature levels during vacuum operation 
produced failures m the metal-to-ceramic seals required at 
the points where the heater coils pass through the quartz outer 
shell of the prototype. The area of this problem is designated 
as area "C" m Fig. 34 . 

B) Preprototype Thruster 

As a result of these problems, a second laboratory thruster was 
fabricated utilizing a metallic external (molybdenum) shell and with 
elimination of as many seals as possible. In fabricating this device 
which was denoted as a "preprototype thruster, " insulation was placed 
between the thruster chamber and the external shell, a configuration "D" 
injector system was incorporated to move the reaction zone away from the 
thruster base plate, and the device was made sufficiently rugged so as to 
make it possible for the unit to be disassembled and reassembled without 
fear of damage A drawing of this device is presented m Fig 35 The 
stainless steel injection tube employed m this device was a commercially 
available hypodermic tube (0.025 in o.d. x 0 015 m. 1. d. ) having an in- 
ternal diameter less than the quenching distance of hydrazine (0.045 m ) 
as reported by Gerstein. 

Eicperimental performance data (characteristic velocity) were 
initially obtained with the preprototype thruster over the range of 
1 x 10 _ 5 to 5 x 10~5 Ib/sec propellant flow rates to provide a direct 
comparison between this unit and the laboratory devices utilized pre- 
viously These results are illustrated m Figs. 36 and 37 with the per- 
formance shown as characteristic velocity, C*, as a function of power 
addition at a propellant flow rate of 5 x 10“^ lb/ sec It is immediately 

obvious that the new prototype gives a considerably better performance 
than the previous prototype. For example, m the unpowered mode of 
operation, the new thruster had a char acteristis tic velocity of 3400 ft/sec 
as compared to 2900 ft/ sec with the older device. It is also interesting 
to note that the characteristic velocity during unpowered operation with 
the new thruster is greater than that obtained with the old device even 
with the addition of B0 watts of power. During unpowered operation, the 
new prototype would have a specific impulse of 172 sec and a delivered 
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3 

thrust of 8 7 x 10" lb^ with an assumed thrust coefficient of 1 65 Hence, 
a thrust level of 10 millipounds could be obtained by increasing the pro- 
pellant flow to a value of 5 71 x 10"^ lb m /sec 

During powered operation, a power addition of 41 watts resulted in 
an increase in the characteristic velocity from 3400 ft/sec to a value of 
3750 ft/sec, and the specific impulse to a value of 192 sec from 172 sec 
obtained during unpowered operation 


Additional data were obtained during unpowered operation at flow 
rates over the range of 1 5 x 10“® to 10 x 10“5 lb/sec as shown in Fig 38 
m the form of characteristic velocity as a function of propellant flow rate 
The characteristic velocity was found to increase monotomcally over the 
entire flow range investigated varying from approximately 2000 ft/sec 
to 4000 ft/sec at 10 x 10 lb/sec at feed rates of 1 5 x 10“^ and 10 x 10 - ^ 
Ib/sec, respectively This behavior is characteristic of a system with 
heat losses which are dependent on the propellant flow rate as can be 
demonstrated by the following analysis. A material balance on the thruster 
may be written in terms of molar flow rates as 


- n - r„V 


(53) 


where r m is the molar generation rate per unit volume (V)* For the 
chemical reaction 

R — "l)P (54) 

The molar generation is related to the reaction rate of the reactant* 

(55) 

For a first order reaction 

r= KLC r 1 = k LC*1 ®<p(-e/rt) , (56) 

with the concentration of the reactants given by 

C R = ' ( 57 ) 

RT 

with the mole fraction (XR) is defined as 

X = - l ~^ - ! (58) 

R v>.+*v i-s+^y 

where the quantity y is the fraction of reactants converted to products 
The input and output molar flow rates are also related for any fraction 
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( 60 ) 


of propellant converted to products since 

r>„ = (i-y+ 

Therefore, 

y ro-0 ft* * (u-o rV 


y n t s rV 


where 


i-a 1 _p 

k exp(-E/l 

L ‘-9 +T>yjRT 

1 V 


An energy balance on the thruster may be written in terms of mass flow 
rates as 

ir> my j-AH r J+ KAy|H w -Hj = (64) 

where the term, — w ( represents the energy losses or addition to 

the system. Letting 

^ — H 0 




-AWr) 


'J* 

» 

or 


J 

Cp 1 

j 


y - • 

[h-Ho' 

}_ 

' H*-U] 



l~AH r 

J mCp 1 

. -aI-WJ 


Substituting equation (62) into the above expression yields 


a-Uc 

-aW*. 


T — * f H w - V) 

J mCp l -AWr J 


where 


Yl *' a m / Ml (68) 

with Mj being the molecular weight of the inlet flow to the system Therefore 

rVI-V _ [ H-H. 1 . WAw f jW-H ] , 6 o, 

m t -A Hr i W>Ga L -AWr J 
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which upon substitution of equation (63) for the reaction rate, r 
yields 


PVM< 1 

( i-v 1 

. k exp| 


iU-U,'] 

l _ hAw ] 

i Hw-Ul 

■RTm 1 


UtJ j 

L-aUv. j 

[ mCp 

L -aHr J 


( 70 ) 


It is also noted that the mass flow through the system can be related to 
the pressure and nozzle exit area since 


w = PA H 


si 


it' 


9*M/ z 

RT \*h J 


where the molecular weight of the exit gases is 

M - Kv 

which implies 

PVH, _ V ( 

Am 


'CART 


-Z V£j R~T ~ 
*•« ) t^T 

A speed of sound Q. v may be defined as 




Hence, equation (70) may be written as 

V kexp(-g/rer) 


U-H, 


Am ^ f 2. ~ ]*Fo -a H r 

or by use of equation (66) 

V k exp(-£/RT) y ( i-y-v-TjyV' 3 ' 


Au 
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wCc 


r u w -u 
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(71) 


(72) 


(73) 


(74) 


(75) 


(76) 


which relates the fraction of the inlet flow converted to reaction products 
and the characteristic length of the thruster since 


L?s V 
A h 


(77) 


The effects of heat losses from the system (Hw H) can be best seen from 
equation (75) 


L* k**p(-EAnO 1-9 _ U~Uo + WA* ( u-wj 

CL l f -AUv VyiC^, (. -A. Hr , 


(78) 
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which may be written as 

:' - & hA L»Kexp(-g/W) (l-U) 

U-H^J (^j)^ (i-^-oy)* 4 - 


\ + 



‘(79) 


or 


F-f (tV ~ 1 i^a _ 

(i-vjvuy)'*. 


\+ / H- VA W \ 

\V\Cp \ H- VAo / 


(80) 


As the losses increase, the characteristic temperature of the system 
decreases and f(T) becomes smaller Hence the quantity 

„ (8! 

/i-a+Tjy)^ 

must increase. Since f(y) decreases momtomcally with increasing values 
of the fraction of propellant decomposed, increasing heat losses implies 
a lower conversion of the feed to reaction pro^ducts. 



Equation(80) may also be employed to demonstrate that in a system 
which experiences heat losses, increasing the flow through the system can 
result m a greater amount of propellant decomposition thereby producing 
an increase m the characteristic temperature of the system. For simplicity, 
consider the system to have a constant wall temperature The non-dimensional 
heat loss term m equation (80) may be written as 

h Aw ( ^ ~ Wwl — tza ^ =? Cy* ( 82 ) 

Vyi Cp ( V-A - vx. ( V\ - P \ jp 


where C n is the Stanton number and A-f is the cross-sectional area for fluid 
flow through the system By definition 



Klq 

Re Pr 


(83) 


with 

Nu= K Re'" Pr™ 


(84) 


with the exponents n and m being less than unity The heat losses may 
then be expressed as 


hAiw ( H-Hw) 
™Gp ( H - H 0 ) 


Avi { 

Af 


K 




(84) 
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Since the Reynolds number is directly proportional to mass flow rate 
through the system 


Re = Patk - W\De 

yAX 

For laminar flow through the system (i.e, n=0 
may be written as 


[ftT)1 llj - 14- KAw 

('~3 + ^3) ,/5C A* 


( 86 ) 


50, m=0 33), equation (80) 



(87) 


This implies that increasing the mass flow decreases the dimensionless 
heat losses and results in an increase m characteristic temperature and 
the fraction of propellant converted to reaction products 


The characteristic velocity data were utilized to calculate specific 
impulse and thrust using an assumed value of 1 65 for the thrust coefficient 
The results of the computations are shown m Fig. 39 Specific impulse 
is found to vary from 120 sec at 2 x 10"® lb m /sec flow rate to 210 sec at 
8 x 10"® lb m /sec. The thrust level developed by the thruster varied from 
2 4 x 10"3 ibjf to 16 8 x 10“-® lbf over the same flow range This is also a 
considerable increase m performance over that obtained with the con- 
figuration "C" thruster as shown m Tables I and II 

Data obtained by Schreib and Pugmire(®) with a configuration "C" 
thruster are included m Table I, whereas, the performance listed m Table II 
for this type of device is a composite of the two sets of data 


A similar series of experiments were conducted to investigate 
thruster performance with power addition to the thruster heater coil using 
propellant flow rates of 2 x 10 to 5 x 10 - 8 lbm/sec, and power inputs of 
20 to 80 watts. The results of these tests are illustrated m Figs. 40-42 . 
At a flow rate of 4 x 10 - ^ lbm/sec, power addition amounting to 80 watts 
was found to produce a 33 percent increase m thruster performance This 
is equivalent to a power requirement of 34 watts for each millipound of 
thrust enhancement Comparison of this result with predicted thrust en- 
hancement curves indicates that the power addition process vanes from 
43 percent efficiency if constant chemical composition is assumed to be 55 
percent efficient if the composition of the reaction products is considered 
variable 


The heater current and voltage measurements obtained during the 
experimental determination of the effects of power addition on thruster 
performance can be employed to provide some indication as to the 
average heater coil temperature Since the heater coil is cemented to 
the internal surface of the alumina cylinder which forms the thrust chamber 
walls, this average temperature may be a good representative value for 
the thruster. The average temperature of the thrust chamber wall can be 
calculated from the resistance of the heater coil since the resistance of the 
molybdenum heater coil is linear m temperature The average heater coil 
temperature was then employed to calculate a characteristic velocity of 
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Table I 


Comparison of Thruster Performance 


Hydrazine Feed 

' 

Characteristic Velocity, ft/ sec 

Rate Ibm/sec 

Configuration 

,r c’' Thruster 

Pre -prototype 


Ref 7 

Thruster 

2 x 10" 5 

1900 


2350 

3 x 10" 5 

2300 


3025 

4 x 10" 5 

2650 


3450 

5 x 10" 5 

2900 

2800 

3700 

6 x 10~5 


2900 

3900 

7 x 10" 5 


3000 

4000 

8 x 10" 5 


3100 

4100 


Table II 

Comparison of Thruster Performance 


Hydrazine 

Configuration 11 c M Thruster 

Pre-protptype Thruster 

Flow Rate 

Isp 

Thrust 

Isp 

Thrust 

Ibm/sec 

sec 

lbf 

sec 

lb f 

* 

2 x 10- 5 

97 

1 9 x 10~ 3 

120 

2 4 x lO" 3 

3 x 10 -5 

118 

3 5 x 10" 3 

155 

4 6 x 10-3 

4 x 10- 5 

136 

5 4 x 10“ 3 

177 

7 1 x lO" 3 

5 x 10" 5 

143 

7.2 x 10" 3 

190 

9.5 x 10- 3 

6 x 10 -5 

149 

8.9 x lO -3 

200 

12 0 x 10" 3 

7 x 10' 5 

154 

10 8 x 10** 3 

206 

i 

14.4 x 10“ 3 

8 x 10~ 5 

159 

12 7 x 10" 3 

210 

| 

16. 8 x 10“ 3 
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the gas decomposition products The results are illustrated m Fig. 43 
which gives a comparison of experimental chaiactenstic velocities 


C* — q PA m 

m 


and calculated values obtained from the theoretical relationship 

qiRT 





TFT 

t-1 


(89) 


where T is the average thrust chamber wall temperature obtained from 
resistance measurements and the specific heat ratio is that for hydrazine 
decomposition products with 40 percent ammonia dissociation 


A somewhat similar procedure was followed to establish average 
thrust chamber wall temperature m conducting those experiments on 
which the results listed in Table II of this report were based However, 
m these tests the voltage to the heater coil was maintained below a level 
of one volt with the power dissipated in the coil being less than one watt 
On the basis of the effects of power addition illustrated in Fig 7 , it is 
evident that this low power dissipation should have little or no effect 
on characteristic velocity The current and voltage to the coil at various 
mass flow rates was processed m a manner similar to that discussed 
above to obtain average thrust chamber wall temperatures and then 
characteristic velocities which could be compared with experimental data 
(Fig 44). 


The calculated characteristic velocity was found to be m good 
agreement with the experimental data m the tests which were conducted 
m what is essentially the unpowered mode of operation (Fig. 44). The 
maximum deviation between the calculated and experimental characteristic 
velocities amounts to approximately 250 ft/sec at propellant flow rates 
m the vicinity of 8 to 10 x 10“^ lbm/sec, where the experimental results 
are nearly constant at a value of 3920 ft/sec. The calculated values are 
m much poorer agreement with experiment in. the powered mode of 
operation (Fig. 43 ) an d the deviation was found to vary quite strongly 
with propellant flow rate Closer examination of ther esults indicate 
that, at the power addition levels considered, the average thrust chamber 
wall temperature is only slightly dependent on flow rate, i e , the wall 
temperature appears to be just a function of the power addition level 


It is also of interest to consider how much the actual temperature 
distribution along the interior surface of the thrust chamber wall differs 
from the mean value For the purposes of this discussion, the thrust 
chamber wall is considered to be an exterior insulated homogeneous 
hollow cylinder which is subjected to a uniform heat flux on its interior 
surface Heat losses from the body were assumed to occur at the surfaces 
at the two ends of the body The steady state temperature distribution m 
such a body can be obtained from the solution of the partial differential 
equation 
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<£l + 1 SE + £l =o 

ir* r hr 


subject to the boundary conditions 


H=o 

r='r 0 



-KSI - q 




U - O 

Z= T 

*Z v 


-kST . 

z=o,u 
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(91) 

(92) 

(93) 


where the boundary condition at the mid-plane of the body reflects the 
requirement that the temperature distribution by symmetrical The 
solution to these equations is simply ' ' 


T(r,lK(r„o) 






(94) 


where the temperature, T (r Q , 0), is the lowest temperature within 
the body The temperature along the interior surface of the hollow 
cylinder is then 


Kn>c>) 


* I * 

zK(«i -0 



where 


(95) 


k<o,q)- 



(96) 


The average temperature under uniform heating conditions may be 
obtained by integrating equation ( 96 ) to give 

.1 

(97) 

n v’ i i >»/ / v i . / i i \ 

T 




T(r„,-z)- l (r w ,o) 


6K(r;-Tf; 


( 98 ) 
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The maximum temperature within the body occurs at the mid-plane 

(Z=Lr/2) 

~T[r Lj z) 

MAX 

Thereby indicating that 

T (tyz) ~T(r,,p) 

or 

T(r w z)= 


= vjb 

4K (n*-r*) 
~ O 667 


( 99 ) 


( 100 ) 


( 101 ) 


The maximum temperature which can occur within the thrust chamber walls 
during unpowered operation is the adiabatic flame temperature of hydrazine 
decomposition products. Assuming 40 percent ammonia dissociation m the 
reaction products, this indicates that the maximum value of "7* ( r v L/g) 
is 2i00°R. A reasonable assumption as to the temperature at 
the end of the thrust chamber walls is 700°R which would imply that the 
average temperature could vary between 67 percent and 78 percent of the 
maximum temperature of the tube wall 

The exact relationship between average wall temperature and exhaust 
gas temperature has not been established. However, the experimental 
data indicates that the average wall temper atur e gives a reasonable estimate 
of the temperature of hydrazine decomposition products when the thruster 
is operated m the unpowered mode. 

In the higher portions of the propellant feed range investigated, flow 
instabilities were observed Typical experimental data obtained which 
demonstrates this phenomenon are presented m Fig. 45 for the thruster 
operated m the unpowered mode. Notice that as the flow rate is increased to 
value of 3 x 10“ 5 Ibm/sec, short duration pulses are observed in the trans- 
ducer which monitors the pressure m the propellant feed line immediately 
upstream of the thruster. No fluctuations are observed, however, m the 
propellant flow rate or thruster chamber pressure. As the hydrazine feed 
is increased, however, the inlet line pressure fluctuations become con- 
tinuous m nature and at high flow rates somewhat random variations in the 
mean value of the inlet pressure are also observed. At the same time, 
increasing the mean flow rate results m fluctuations in propellant flow and 
m the chamber pressure. The maximum pressure fluctuations observed 
m these experiments amount to approximately +3 percent of the mean 
nozzle pressure while the mass flow rate fluctuations were on the order of 
+5 percent 
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These instabilities can be explained as follows Consider 
the porous zircoma and feed tube to be a homogeneous body which 
has a uniformly distributed volumetric heat release rate (Q) m that 
portion of body downstream of the end of feed tube as shown schematically 
m Fig. 46. Furthermore, let us assume that there are no losses from 
the sides of the body or from the downstream end while the upstream end 
of the tube is maintained at a constant temperature, T . The temperature 
distribution m that portion of the body downstream of the feed tube, 1 . e. 
region 1, in Fig. 46, is given as (13) 


eo 


T-X = 



cosj" ta* p ^ 'TV 


YVsO 



( 102 ) 


which at steady state yields 


T-T b 


= OL 
aK 



(103) 


The heat which is transferred from the reaction zone to the remaining 
portion of the body is simply 




= Q.L 


I 


(104) 


This energy is absorbed by the propellant m the upstream area of the 
body m the process of vaporization or 



It is evident that at some flow rate the energy transferred down the tube 
is less than that required to vaporize the propellant since the maximum 
temperature which can be obtained in reaction region is the adiabatic flame 
temperature of the decomposition reaction. The maximum flow for which 
all the fluid is vaporized is then 




|c p (T b -X)+£>U 





2-k ( 'a? 


(106) 


or 




ZK(Tqg-Tb^ 

L, |C P K- 


(107) 
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At flow rates m excess of m^, the fluid will be only partially vaporized 
before it enters the reaction zone, while at lower flow rates, the pro- 
pellant will be completely vaporized. Employing a thermal conductivity 
of 0 004 cal/cm-K-sec and the adiabatic flame temperature corresponding 
to an ammonia dissociation fraction of 0 40, the flow rate calculated from 
Equation (107) is 2.0 x 10“ 5 lb m /sec 

At low flow rates, we would expect the propellant flow to be 
completely vaporized m the upstream portion of the tube and hence, 
would have the characteristics of a gaseous feed system while at high 
flows it would behave as a liquid injection system This agrees well 
with what is observed experimentally as is shown by the following dis- 
cussion At low flows, the propellant feed is directly proportional to 
the inlet pressure immediately downstream of the propellant control valve 
as shown in Fig 47, whereas at higher flow rates, the mass flow is non- 
linear m the inlet pressure These data were also correlated with an ex- 
pression of the form 

m= PP, (108) 


where the quantity F depends upon the ratio of inlet pressure to 
thruster chamber pressure 



for *=1 30 


(109) 


-<056 (110) 


These expressions are essentially the relationships which define the 
compressible gas flow through an orifice The correlation of the mass 
flow withethe inlet pressure is presented in Fig 48. At low flows, it 
is seen that 


Yn= KFP, 


\ 

\ 


(ill) 


as would be expected with a gaseous feed system while at high flow rates, 
the deviation m the observed flow from that obtained by extrapolation of 
the linear portion of the curve increases with increasing flow rate This 
deviation is believed to be the result of the propellant entering the reaction 
zone in a partially vaporized condition The instabilities inherent in a 
two -phase flow system are well documented^* and it is believed that 
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this is the cause of the pressure fluctuations observed in the preprototype 
thruster. 


There is additional data which also supports this conclusion During 
operation of the preprototype thruster m the powered mode, it was observed 
that at any operating point where there were instabilities, increasing the 
power input to the heater coil decreases the magnitude of the fluctuations. 
Typical recorder traces which demonstrate this effect are shown in Fig 49 
Obviously, some portion of the power applied to the coil is transferred 
to the zircoma and in turn results m an increase in the propellant flow rate 
at which two-phase flow first occurs 


A second series of experiments were conducted with the preprototype 
thruster on a thrust and impulse measuring system The thrust and impulse 
measuring system consists of a "frictionless " inertial table whose angular 
velocity is measured by a rate gyro The table is suspended on a multi- 
filament steel cable. The wire is mounted to the inner race of a thrust 
bearing A servo-torque motor slaves the bearing position to the table 
position as a function of the error signal generated by the misalignment of 
a polarized light source (on the table) and sensor (on the bearing mount) 
Thus, the torque due to the angular deflection of the support wire may 
be maintained at a value less than that equivalent to 1 x 10“^ lb force act- 
ing at a two foot radius. An artist's representation of the thrust measur- 
ing system and its associated vacuum system are illustrated m Fig. 50 

The application of Newton's second law to the table results in the 
following equation 

v* 

^ ( Fr) h = ( 112 ) 

where 

F= external tangential force (impulse thrust) 

R=radius of tangential component of F 
J=moment of inertia of the table 

©=rate of change of angular velocity 


This measurement system may accommodate any n um ber of torques, 
generally, a single torque is considered However, couples, as a pair of 
yaw thrusters, may be studied individualy or as a couple. 

For impulse measurements of a single thruster, equation (112) is 
integrated once , 


r 

Fits 1= j(0 t -0> 


(113) 
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Fxg 49 Effect of Power Addition on Thruster Instabilities (i$i=5 0x10 “^lb m / sec) 
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The gyro output is described by 


© = K ( E wv ) {114) 

where 

E mv = millivolt output signal 
k = proportionality constants 

Substitution of equation (114) into equation (113) results m the working 
equation for impulse measurements 



The constant (Jk) is determined by calibration tests. The table 
is calibrated by applying known torques to the table and monitoring the 
gyro output during a calibration test. The applied force m equation (115) 
is constant so F may be brought outside the integral sign. Equation (116) 
describes the calibration test 

r c - JK a (Hk) 

At 

where subscript c refers to calibration parameters 


Rearrangement of equation (116) results m the working equation for 
evaluating (Jk) of equation (115) 

C s JK = t\ — ( 117 ) 

\ A /c 

The total impulse delivered per flow valve pulse as calculated from 
equation (115) includes all losses due to pressure dependent, specific im- 
pulse, and nozzle viscous losses. If sufficient pulses are made, the total 
variation in impulse bits will include the variation m the plenum pressure 
within the deadband limits of the regulator. 


It is noteworthy to point out some experimentally significant 
characteristics of equations (113) and{117). All the parameters of equation 
(117) are readily measurable and no compounding of inaccuracies result 
from its application to the calibration testing. The time different (At), 
and the millivolt difference ( Amv), may be determined to within 5% 
The mass of the test weight and the moment arm may be determined to 
even greater accuracies. Thus, the evaluation of K should be within 
1%. The application of equation (115) to the impulse mode of testing re- 
sults m a direct measurement of impulse delivered as determined by a 


79 



voltage difference No further data analysis (as integration of force- 
time responses) are required. Thus, good resolution of the desired 
parameter, impulse bit, is obtained. Due to the nature of the rate 
gyro, frequency response is not a limiting factor The impulse 
measurement is dependent on resolving the difference in angular rate 
before and after the application of the impulsive thrust, and not on 
the ability of the gyro to faithfully tract the table velocity during the 
application of the impulsive force Thus, good integrated impulse 
resolution is obtained at the sacrifice of the detailed thrust-time history. 


A schematic representation of a typical gyro output trace is 
shown m Fig. 51 . Chronologically, the following sequence of events 
occur during a typical cycle of the gyro output signal At time t 0 , 
the thrust delivered by the thruster is causing the table to accelerate 
in a counter-clockwise direction although the table is at rest. At time 
tp the table is moving at its maximum velocity and an auxiliary boom 
on the thruster engages an elastic bumper which changes the motion 
of the table from the counter-clockwise to the clockwise direction. 

This is accomplished during the period tj to tz During the time be- 
tween t 2 and 13 , the table is being decelerated at a constant rate by the 
thruster The table comes to rest at time t 3 and is accelerated towards 
the bumper until time t 4 Referring to equation (115), the thrust is 
determined by the slope of the gyro output between times t 2 and t 4 . It 
is noted, however, that this representation is idealized m that the gyro 
output is strictly a linear function of time In actual experiments, how- 
ever, the output is not exactly linear as is shown by the experimental 
data illustrated m Fig. 52. These anomalies in the slope are caused 
by extraneous phenomena within the vacuum tank such as recirculation 
flow patterns, etc. Since most effects tend to produce a thrust decrement, 
it can be argued that the most logical slope of the gyro output curve to use 
in calculation of thrust is the maximum slope. For the example illustrated 
m Fig. 52, this- amounts to a thrust obtained from the maximum slope which 
is 9 percent higher than that would be obtained from the average slope. 


Steady state data obtained with the preprototype thruster are 
illustrated in Fig. 53 . Data were reduced using both the maximum and 
average slope of the gyro output curve. It is seen that the deviation 
between the two methods amounts to approximately 9 percent. In order 
to provide a conservative evaluation of thruster performance, the average 
slope method has been adopted for representation of the experimental data 

Characteristic velocity and thrust coefficients obtained from ex- 
perimental data are depicted m Fig. 54. Characteristic velocity is m 
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close agreement with similar data discussed previously Thrust 
coefficients were determined from the definition 



where the thrust employed was that obtained from the average slope of 
wire table gyro output curve The thrust coefficient was found to remain 
relatively constant with a mean value of 1 55 over the entire range of 
propellant flow rates investigated If the maximum slope of the gyro 
output curve were used to determine thrust, the average thrust coefficient 
would be increased to a value of 1 70. 

A similar series of experiments were performed with a second 
preprototype thruster which had the porous insulator indicated m Fig 35 
replaced with a solid refractory insulator (lava) such that the characteristic 
length of the thrust chamber was identical for the two devices. During 
operation of this unit on the thrust table, improvements were made to the 
propellant feed system which made it possible to increase the propellant 
supply pressure from the previous maximum value of 80 psid to 180 psid 
Hence, thrust data could be obtained over a wider range of propellant flow 
rates (3 x 10" 5 to 25 x 10 -5 lb/sec). Steady state thruster performance 
obtained with this unit is illustrated in Figs 55,56 Thrust was found to vary 
from 5 millipounds at the minimum flow rate to 59 millipounds at the high- 
est propellant flow investigated using the average slope of gyro output 
signal The corresponding specific impulse increased momtonically 
from a value of 160 sec. to 238 sec when propellant feed from 3 x 10"^ lb/ sec 
Comparison of the performance of the two preprototype thrusters indicate s 
that, as expected, the porous insulator was more effective in reducing energy 
losses from the thruster since the lava insulator has the higher thermal con- 
ductivity 

Pulsed mode operational data were obtained with the preprototype 
thruster containing the lava insulator for duty cycles varying from 1 to 50 
percent duty cycle with a 1 sec thrust period These investigations 
were performed with the unit operating m both the unpowered and powered 
modes The experiments were performed with this unit since it was ex- 
pected that it would have a response time which would be considerably shorter 
than that of the preprototype containing the porous insulator. In the original 
preprototype, steady state nozzle pressure is achieved only after both the 
thrust chamber and the volume occupied by the porous insulator reach their 
steady state values, whereas, with the solid insulator only the thrust chamber 
needs to be brought up to pressure to attain steady state 

Pulse data were gathered at nominal steady state thrust levels of 
9 x 10"3 Ibf and 17 5 x 10 Ibf. In performing these experiments, the flow 
regulating valve was adjusted to give the steady state flow rate required to 
obtain the desired thrust level The flow regulating valve setting was main- 
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FIG. 56 MODIFIED PRE PROTOTYPE THRUSTER PERFORMANCE , POWERED 

OPERATION 





tamed m this position during the pulsed mode of operation Typical 
recorder traces obtained m these experiments are illustrated m 
Figs. 57-39 Upon examination of these data several phenomena 
become apparent, namely. 

a) There is a considerable lag in the output of the mass flow 
measuring system pressure transducer as demonstrated by the long 
period required for mass flow signal to return to zero after the pro- 
pellant valve is shut 

b) The time required for the thrust to decay to zero is on the 
order of 600 sec. 


c) The maximum propellant flow rate during any pulse is in 
excess of that obtained during steady state operation. 


Therefore, specific impulse obtained at any duty cycle could be 
obtained only after integration of the mass flow rate signal 

T = j F dt 

-*-sp 


(119) 


where 


jFdt = IK(E - E mVi ) 


( 120 ) 


Pulse mode operation data are illustrated m Figs. 60-62 for both 
powered and unpowered modes of operation m normalized with respect to 
parameters obtained at a 50 percent duty cycle. It is noted that at each 
thrust level the specific impulse obtained decreased quite rapidly with 
decreasing percent duty cycle to a value of 40-45 percent of the value ob- 
tained at a 50 percent duty cycle This was accompanied by an increasing 
integrated propellant flow rate and impulse bit as percent duty cycle was 
decreased from 50 percent to 1 percent. The observed drop m specific 
impulse is considerably greater than similar data obtained with catalytic 
hydrazine thrusters (16) For example, it is reported that with a 
no mi nal 5 pound thruster a degradation of specific impulse to 78 percent 
and 71 percent of steady state were observed when operating at conditions 
corresponding to 4 2 lbf and 2 5 lb| with a 0 4 percent duty cycle 


A limited amount of data were obtained with the preprototype 
thruster ope rating with a pulse width of less than 1 sec. These results 
are listed m Table III, below 
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PULSE MODE RECORDER TRACE'S , MODIFIED 
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FIG. 60 MODIFIED PRE PROTOTYPE THRUSTER PERFORMANCE AT 

50 % DUTY CYCLE ( I SECOND PULSE ) 





Ft<3. 61 EFFECT OP PUTY CYCLE ON MODIF1EO PRE PROTOT Y PET THPUSTE R 

PERFORMANCE 








FIG <b2. EPFECT OF "DUTY CYCLE ON MODIFIED PREPROTOTYPE 

THRUSTER PERFORMANCE. 
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Table III 

Summary of Short Pulse Data 
Preprototype Thruster, Lava Insulator 


Pulse Width 
sec. 

% Duty 
Cycle 

Specific Impulse 
sec 

Impulse Bit 
(lbf - sec) x 10^ 

O 0 

100 

228 

49 x 10" 2 3 

0 10 

10 

200 

21 x 10" 3 

0.03 

3 

189 

15 x 10~ 3 

CO 

100 

238 

59 x 10' 3 

0 10 

1 

180 

41 x 10" 3 


Although these data do not provide sufficient definition of performance with 
short pulses, the degradation in specific impulse is considerably less 
(I S p =o 75 l S p ss ) at a 1 percent duty cycle than that observed with the 1 sec. 
pulse. 


It is believed that the reason for the discrepancy in the two sets of 
data can be explained by means of the following hypothesis as to the sequence 
of events occurring within the thruster after opening the propellant feed valve 

1) At the time at which the propellant feed valve is opened, the 
propellant feed tube within the thruster has a temperature greater than 

the boiling point of the liquid propellant and the pressure within the thruster 
is low As propellant flows down the feed tube, it is partially vaporized and 
the propellant flow rate is relatively low (when compared to that obtained 
later m the operating period) since the pressure drop per unit mass flow is 
high for a two phase flow During this period, the temperature of the inlet 
tube is continuously decreasing until it drops to a value below the boiling 
point of the propellant In this part of the operating cycle, the specific im- 
pulse obtained is hypothesized to be close to that which would he achieved m 
steady state operation 

2 ) During the second portion of the thrust period, propellant enters 
the thruster in the liquid state For given values of supply and thrust chamber 
pressure, the propellant feed rate achieved for liquid flow is considerably 

hi gher than with a two phase mixture. Hence, there is a rapid increase in 
propellant flow to values greater than that achieved m either the first portion 
of the pulse or during steady state operation The net effect of this high flow 
rate is to partially flood the thruster which results in thrust being delivered 
at specific impulses characteristic of partially reacted hydrazine These 
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hypotheses are supported by experimental pressure traces obtained 
for long pulses (see Figs 26, 27 and 28) which indicate that steady 
state behavior would be achieved only with pulses much longer than 
1 sec duration. 

3) Following shutting of the propellant feed valve at the end of 
the pulse, thrust and thrust chamber pressure decrease to a low value. 
However, thrust continues to be delivered as a fluid contained within 
the ullage volume between the valve seat and the thrust chamber vaporizes. 
During this period, the temperature of the propellant feed tube increases to 
attain values m excess of the propellant boiling point until the time at which 
the flow valve is opened at the start of the next pulse. 

This hypothetical sequence of events would explain the difference 
between the observed behavior of the thruster when operated with short 
pulses (less than 0 10 sec.) and long pulses (1.0 sec.) However, it 
must remain strictly conjecture since the experimental data cannot be 
reduced to illustrate thrust as a function of time due to the lag m the mass 
flow measuring system signal and the inability of the thrust measuring 
system to accurately follow the variation of thrust with time 

V Materials Compatibility 

In the course of these investigations, several materials problem 
areas came to light, namely 

a) The effects of hydrazine vapor and/or decomposition products 
on stainless steel at high temperatures 

b) The effects of liquid hydrazine storage m stainless steel 

lines. 

c) The anomalous behavior of stabilized zircoma when employed 
as a flow distributor m the hydrazine resistojet. 

These problem areas are discussed below. 

(1) During operation of the preprototype thruster m the pulsed 
mode, it was observed that it was impossible to remove the stainless steel 
feed line from the porous zircoma tube. The thruster was disassembled 
and the injector assembly was sectioned m the manner shown schematically 
m Fig 63 , and photo micrographs were taken of the various cross -sections 
(Figs. 64-66). In the two photo micrographs, which include both the ceramic 
and the metal tube (Fig. 61), degradation of the tube by hydrazine is readily 
apparent with the amount of corrosion increasing with distance down the tube. 
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At the downstream end of the tube, there is corrosion of some degree 
throughout the entire tube wall, this degradation was accompanied by 
a decrease in the inner diameter and an increase m the outer diameter 
of the tube This deformation has occurred to such an extent that the 
outer surface of the feed line has deformed m those regions where it is 
m contact with the sourrounding porous ceramic. Downstream of the end - 
of the feed tube the porous zirconia appears to be physically similar to 
virgin material {see Fig. 66) with the exception of small pieces of metal 
m the passage down the axis of the ceramic It is not known whether these 
particles were a direct result of chemical attack on the stainless steel or 
were produced in the process of cross -sectioning the assembly Since 
the particular unit undeF'discuSsion was used for a total run time in the 
vicinity of 10 hours, it is apparent that corrosion of stainless steel by 
hydrazine can be a serious problem 

The mounted specimens shown in these photo micrographs and 
a sample of the stainless steel feed line have been forwarded to the 
NASA Goddard Space Flight Center for evaluation On the basis of 
metallographic examination and microhardness test results, it 
was concluded that both the inside and outside surfaces of the stainless 
steel feed tube were mtrided during operation of the thruster On the 
basis of these investigations, it was recommended that Inconel or tantalum 
should be employed m fabrication of the feed tube or that the feed tube 
be plated to minimize nitriding with gold being the most likely candidate 

(2) Operation of the liquid hydrazine feed system in these in- 
vestigations disclosed that there was a distinct possibility that liquid 
hydrazine is slowly decomposed on steel surfaces Evidence as to this 
phenomenon is as follows It was observed throughout the program that 
the response of the flow measuring was far removed from that which would 
be expected with a pure liquid system as demonstrated by lags m transient 
response of the differential pressure transducer across the flow metering 
orifice The slow transient response of the system can only be explained 
by compliance within the feed lines m the form of gas pockets. Direct 
experimental evidence of gas generation within the system was obtained 
from the necessity of periodic venting of the gas lines m the system How- 
ever, sifter venting the lines it was observed that the compliance gradually 
increased with time. 

(3) Anomalous behavior of the stabilized zirconia flow distributor 
was observed throughout the program. During steady state operation of 
both laboratory and preprototype thrusters, the use of stabilized zirconia 
was found to present little or no difficulty except when power addition at 
high levels {60-80 watts) was attempted. At the high power levels, the 
zirconia tube was found to be cracked upon disassembly of the thruster 
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This behavior was attributed to increases m radial temperature gradient 
in the thruster and hence, larger thermal stresses in the powered mode 
of operation. During pulsed mode of operation, the zircoma was found 
to crack with greater frequency probably due to high thermal stresses 
generated by the sudden surges in propellant flow rate Although this is 
the most likely reason for the observed behavior, the possibility of 
leaching of the calcia stabilizing agent by hydrazine should not be eliminat- 
ed as the reason for the cracking. Experimental evaluation of a thruster 
with a cracked zircoma tube indicated that the cracks did not seriously 
effect thruster performance m that it produced only a 4 percent decrease 
m characteristic velocity However, the distinct possibility of ceramic 
particles lodging within the thruster nozzle throat with the accompanying 
decrease in performance due to nozzle inefficiencies leads to the con- 
clusions that a porous metal structured be considered m place of the 
ceramic The most likely candidate materials for this flow distributor 
are Inconel and tantalum Although the presence of metal surfaces with 
the thruster will produce some dissociation of ammonia and hence, will 
result m lower specific impulses, it is felt that the dissociation can be 
minimized by decreasing the residence time within the structure. 
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